A unified first-principles study of Gilbert damping, spin-flip diffusion and resistivity 

in transition metal alloys 



Anton A. Starikov, 1 Paul J. Kelly, 1 Arne Brataas, 2 Yaroslav Tserkovnyak, 3 and Gerrit E. W. Bauer 4 

1 Faculty of Science and Technology and MESA + Institute for Nanotechnology, 
University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands 
2 Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 
3 Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA 
4 Kavli Institute of NanoScience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands 

(Dated: May 20, 2011) 

Using a formulation of first-principles scattering theory that includes disorder and spin-orbit 
coupling on an equal footing, we calculate the resistivity p, spin flip diffusion length l B f and the 
Gilbert damping parameter a for Nii-^Fe^ substitutional alloys as a function of x. For the tech- 
nologically important NisoFe2o alloy, permalloy, we calculate values of p = 3.5 ± 0.15 ^Ohm-cm, 
l a f = 5.5 ± 0.3 nm, and a — 0.0046 ± 0.0001 compared to experimental low-temperature values in 
the range 4.2 — 4.8 /iOhm-cm for p, 5.0 — 6.0 nm for l s f, and 0.004 — 0.013 for a indicating that the 
theoretical formalism captures the most important contributions to these parameters. 

PACS numbers: 72.25.Rb, 71.70.Ej, 72.25.Ba, 75.40.Gb, 76.60.Es 



Introduction. The drive to increase the density and 
speed of magnetic forms of data storage has focussed at- 
tention on how magnetization changes in response to ex- 
ternal fields and currents, on shorter length- and time- 
scales pp. The dynamics of a magnetization M in an ef- 
fective magnetic field H e g is commonly described using 
the phenomenological Landau-Lifshitz-Gilbert equation 
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where M s = |M| is the saturation magnetization, G(M) 
is the Gilbert damping parameter (that is in general a 
tensor) and the gyromagnetic ratio 7 = gfis/fi is ex- 
pressed in terms of the Bohr magneton fig and the Lande 
g factor, which is approximately 2 for itinerant ferromag- 
nets. The time decay of a magnetization precession is 
frequently expressed in terms of the dimensionless pa- 
rameter a given by the diagonal element of G/^/M s for 
an isotropic medium. If a non-equilibrium magnetization 
is generated in a disordered metal (for example by inject- 
ing a current through an interface), its spatial decay is 
described by the diffusion equation 
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in terms of the spin accumulation A/i, the difference be- 
tween the spin-dependent electrochemical potentials fi s 
for up and down spins, and the spin-flip diffusion length 
l s { [21 IS]- In spite of the great importance of a and l s { , 
our understanding of the factors that contribute to their 
numerical values is at best sketchy. For clean ferromag- 
netic metals [I] and ordered alloys [5] however, recent 
progress has been made in calculating the Gilbert damp- 
ing using the Torque Correlation Model (TCM) [6] and 
the relaxation time approximation in the framework of 



the Boltzmann equation. Estimating the relaxation time 
for particular materials and scattering mechanisms is in 
general a non-trivial task and application of the TCM to 
non-periodic systems entails many additional complica- 
tions and has not yet been demonstrated. Thus, the the- 
oretical study of Gilbert damping or spin-flip scattering 
in disordered alloys and their calculation for particular 
materials with intrinsic disorder remain open questions. 

Method. In this paper we calculate the resistivity p, 
spin-flip diffusion length ^ s f and Gilbert damping param- 
eter a for substitutional Nii-^Fea; alloys within a single 
first-principles framework. To do so, we have extended a 
scattering formalism [7] based upon the local spin den- 
sity approximation (LSDA) of density functional theory 
(DFT) so that spin-orbit coupling (SOC) and chemical 
disorder are included on an equal footing. Relativistic 
effects are included by using the Pauli Hamiltonian. 

For a disordered region of ferromagnetic (F) alloy sand- 
wiched between leads of non-magnetic (N) material, the 
scattering matrix S relates incoming and outgoing states 
in terms of reflection (r) and transmission matrices (i) 
at the Fermi energy. To calculate the scattering ma- 
trix, we use a "wave-function matching" (WFM) scheme 
7\ implemented with a minimal basis of tight-binding 
linearized muffin-tin orbitals (TB-LMTOs) [8]. Atomic- 
sphere-approximation (ASA) potentials [5] are calculated 
self-consistently using a surface Green's function (SGF) 
method also implemented [9] with TB-LMTOs. Charge 
and spin densities for binary alloy A and B sites are calcu- 
lated using the coherent potential approximation (CPA) 
[TU] generalized to layer structures [5] . For the transmis- 
sion matrix calculation, the resulting spherical potentials 
are assigned randomly to sites in large lateral supercells 
(SC) subject to maintenance of the appropriate concen- 
tration of the alloy [7]. Solving the transport problem 
using lateral supercells makes it possible to go beyond 
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effective medium approximations such as the CPA. Be- 
cause we are interested in the properties of bulk alloys, 
the leads can be chosen for convenience and we use Cu 
leads with a single scattering state for each value of crys- 
tal momentum, ky . The alloy lattice constants are de- 
termined using Vegard's law and the lattice constants of 
the leads are made to match. Though NiFe is fee only 
for the concentration range < x < 0.6, we use the fee 
structure for all values of x. 

For the self-consistent SGF calculations (without 
SOC), the two-dimensional (2D) Brillouin zone (BZ) cor- 
responding to the lxl interface unit cell was sampled 
with a 120 x 120 grid. Transport calculations including 
spin-orbit coupling were performed with a 32 x 32 2D BZ 
grid for a 5 x 5 lateral supercell, which is equivalent to 
a 160 x 160 grid in the 1 x 1 2D BZ. The thickness of 
the ferromagnetic layer ranged from 3 to 200 monolay- 
ers of fee alloy; for the largest thicknesses, the scattering 
region contained more than 5000 atoms. For every thick- 
ness of ferromagnetic alloy, we averaged over a number 
of random disorder configurations; the sample to sample 
spread was small and typically only five configurations 
were necessary. 

Resistivity. We calculate the electrical resistivity to 
illustrate our methodology. In the Landauer-Biittiker 
formalism, the conductance can be expressed in terms of 
the transmission matrix t as G — (e 2 //i)Tr {ttf} [TTlll2j . 
The resistance of the complete system consisting of ideal 
leads sandwiching a layer of ferromagnetic alloy of thick- 
ness L is R(L) = 1/G(L) = l/G S h + 2i? if + R h (L) where 
Gsh = (2e 2 //i) N is the Sharvin conductance of each lead 
with TV conductance channels per spin, i?;f is the interface 
resistance of a single N|F interface, and Rb(L) is the bulk 
resistance of a ferromagnetic layer of thickness L [7J H3] . 
When the ferromagnetic slab is sufficiently thick, Ohmic 
behaviour is recovered whereby Rb{L) ~ pL as shown in 
the inset to Fig. [I] for permalloy (Py = Ni 80 Fe2o) an d 
the bulk resistivity p can be extracted from the slope 
of R(L) For currents parallel and perpendicular to 
the magnetization direction, the resistivities are different 
and have to be calculated separately. The average resis- 
tivity is given by p = (pu + 2pj_)/3, and the anisotropic 
magnetoresistance ratio (AMR) is (py — p±)/p. 

For permalloy we find values of p = 3.5 ± 0.15 ^Ohm- 
cm and AMR = 19 ± 1%, compared to experimental low- 
temperature values in the range 4.2 — 4.8 ^Ohm-cm for 
p and 18% for AMR [15]. The resistivity calculated as a 
function of x is compared to low temperature literature 
values [T3] in Fig. [I] The plateau in the calculated values 
around the Py composition appears to be seen in the 
experiments by Smit and Jaoul et al. [15] . The overall 
agreement with previous calculations is good [16]. In 
spite of the smallness of the SOC, the resistivity of Py 
is underestimated by more than a factor of four when it 
is omitted, underlining its importance for understanding 
transport properties. 
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FIG. 1. Calculated resistivity as a function of the concen- 
tration x for fee Nii-^Fez binary alloys with (solid line) and 
without (dashed-dotted line) SOC. Low temperature experi- 
mental results are shown as symbols [15]. The composition 
Ni8oFe2o is indicated by a vertical dashed line. Inset: resis- 
tance of Cu|Ni8oFe2o|Cu as a function of the thickness of the 
alloy layer. Dots indicate the calculated values averaged over 
five configurations while the solid line is a linear fit. 



Three sources of disorder which have not been taken 
into account here will increase the calculated values of 
p\ short range potential fluctuations that go beyond the 
single site CPA, short range strain fluctuations reflecting 
the differing volumes of Fe and Ni and spin disorder. 
These will be the subject of a later study. 

Gilbert Damping. Recently, Brataas et al. showed 
that the energy loss due to Gilbert damping in an N|F|N 
scattering configuration can be expressed in terms of the 
scattering matrix S [T7]- Using the Landau-Lifshitz- 
Gilbert equation 0, the energy lost by the ferromagnetic 
slab is, 
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where m = M/M, is the unit vector of the magnetization 
direction for the macrospin mode. By equating this en- 
ergy loss to the energy flow into the leads [THj associated 
with "spin-pumping" |19j . 
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the elements of the tensor G can be expressed as 
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Physically, energy is transferred slowly from the spin de- 
grees of freedom to the electronic orbital degrees of free- 
dom from where it is rapidly lost to the phonon degrees 
of freedom. Our calculations focus on the role of elastic 
scattering in the rate-limiting first step. 

Assuming that the Gilbert damping is isotropic for cu- 
bic substitutional alloys and allowing for the enhance- 
ment of the damping due to the F|N interfaces [T5l421j . 
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FIG. 2. Calculated zero temperature (solid line) and exper- 
imental room temperature (symbols) values of the Gilbert 
damping parameter as a function of the concentration x for 
fee Nii-^Fea: binary alloys [2H-23 . Inset: total damping of 
Cu|Ni8oFe2o|Cu as a function of the thickness of the alloy 
layer. Dots indicate the calculated values averaged over five 
configurations while the solid line is a linear fit. 



the total damping in the system with a ferromagnetic slab 
of thickness L can be written G(L) = Gn + Gb{L) where 
we express the bulk damping in terms of the dimension- 
less Gilbert damping parameter Gb(L) = a^M s (L) — 
a"/fx s AL, where p s is the magnetization density and A is 
the cross section. The results of calculations for Ni so Fe2o 
are shown in the inset to Fig. [2] where the derivatives of 
the scattering matrix in ^ were evaluated numerically 
by taking finite differences. The intercept at L = 0, Gif, 
allows us to extract the damping enhancement [20] but 
here we focus on the bulk properties and leave consid- 
eration of the material dependence of the interface en- 
hancement for later study. The value of a determined 
from the slope of G(L)/{jp s A) is 0.0046 ± 0.0001 that 
is at the lower end of the range of values 0.004 — 0.013 
measured at room temperature for Py [2"lTf2"3"] . 

Fig. [2] shows the Gilbert damping parameter as a func- 
tion of x for Nii-^Fe^ binary alloys in the fee structure. 
From a large value for clean Ni, it decreases rapidly to a 
minimum at x ~ 0.65 and then grows again as the limit 
of clean fee Fe is approached. Part of the decrease in 
a with increasing x can be explained by the increase in 
the magnetic moment per atom as we progress from Ni 
to Fe. The large values of a calculated in the dilute al- 
loy limits can be understood in terms of conductivity-like 
enhancement at low temperatures |24) that has been ex- 
plained in terms of intraband scattering [H [B] . The trend 
exhibited by the theoretical a(x) is seen to be reflected 
by experimental room temperature results. In spite of 
a large spread in measured values, these seem to be sys- 
tematically larger than the calculated values. Part of this 
discrepancy can be attributed to an increase in a with 
temperature [2"!?1 |2"B"] . 

Spin diffusioji. When an unpolarized current is in- 
jected from a normal metal into a ferromagnet, the polar- 
ization will return to the value characteristic of the bulk 
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FIG. 3. Fractional spin-current densities for electrons injected 
at z — from Cu into Ni8oFe2o alloy. Calculated values (sym- 
bols) and fits to Eq. (JsJ) (solid lines). 



ferromagnet sufficiently far from the injection point, pro- 
vided there are processes which allow spins to flip. Fol- 
lowing Valet-Fert [3] and assuming there is no spin-flip 
scattering in the N leads, we can express the fractional 
spin current densities p^C-l-) = jt(-f) j J as a function of 
distance z from the interface as 
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where J is the total current through the device, 
and are the currents of majority and minority elec- 
trons, respectively, l s { is the spin-diffusion length, p F — 
(p^ + p^)/4 is the bulk resistivity and /3 is the bulk spin 
asymmetry (p^ — p i )/{p^ + p^). The interface resistance 
rf f = (rf ( + rJ f )/4, the interface resistance asymmetry 

7 = (*"if 4- — r ii)/( r ii + r if) an d the interface spin-relaxation 
expressed through the spin-flip coefficient S [26] must be 
taken into consideration resulting in a finite polarization 
of current injected into the ferromagnet. The correspond- 
ing expressions are plotted as solid lines in Fig. [3] 

To calculate the spin-diffusion length we inject non- 
polarized states from one N lead and probe the transmis- 
sion probability into different spin-channels in the other 
N lead for different thicknesses of the ferromagnet. Fig. [3] 
shows that the calculated values can be fitted using ex- 
pressions Q if we assume that J a /J = G a /G, yielding 
values of the spin-flip diffusion length ^ s f = 5.5 ± 0.3 nm 
and bulk asymmetry parameter j3 — 0.678 ± 0.003 for 
NigoFc2o alloy compared to experimentally estimated val- 
ues of 0.7 ± 0.1 for (3 and in the range 5.0 — 6.0 nm for 

l s { and /3 are shown as a function of concentration x 
in Fig. [4] The convex behaviour of /3 is dominated by 
and tracks the large minority spin resistivity p^ whose 
origin is the large mismatch of the Ni and Fe minority 
spin band structures that leads to a ~ a;(l - x) concen- 
tration dependence of p^ {x) [16]. The majority spin band 
structures match well so p^ is much smaller and changes 
relatively weakly as a function of x. The increase of l s { 
in the clean metal limits corresponds to the increase of 
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FIG. 4. Spin-diffusion length (solid line) and polarization j3 as 
a function of the concentration x for Nii-aFe^ binary alloys. 



the electron momentum and spin-flip scattering times in 
the limit of weak disorder. 

In summary, we have developed a unified DFT-based 
scattering theoretical approach for calculating transport 
parameters of concentrated alloys that depend strongly 
on spin-orbit coupling and disorder and have illustrated 
it with an application to NiFe alloys. Where comparison 
with experiment can be made, the agreement is remark- 
ably good offering the prospect of gaining insight into 
the properties of a host of complex but technologically 
important magnetic materials. 
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